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A critical ingredient for realizing large-scale quantum information processors will be the ability to 
make economical use of qubit control hardware. We demonstrate an extensible strategy for reusing 
control hardware on same-frequency transmon qubits in a circuit QED chip with surface-code¬ 
compatible connectivity. A vector switch matrix enables selective broadcasting of input pulses to 
multiple transmons with individual tailoring of pulse quadratures for each, as required to minimize 
the effects of leakage on weakly anharmonic qubits. Using randomized benchmarking, we compare 
multiple broadcasting strategies that each pass the surface-code error threshold for single-qubit gates. 
In particular, we introduce a selective-broadcasting control strategy using five pulse primitives, which 
allows independent, simultaneous Clifford gates on arbitrary numbers of qubits. 


Building a fault-tolerant quantum computer requires 
the ability to efficiently address and control individual 
qubits in a large-scale system. Many leading exper¬ 
imental quantum information platforms, among them 
trapped ions [T], electronic spins in impurities and quan¬ 
tum dots [2] and superconducting circuits [3], employ 
qubits with level transitions in the microwave frequency 
domain. Addressing these transitions often involves ex¬ 
pensive microwave electronics scaling linearly with the 
number of qubits. To move beyond the state of the art in 
microwave-frequency quantum processors, such as those 
recently used for small-scale quantum error correction in 
superconducting circuits gHS], it will already be bene¬ 
ficial to have a hardware-efficient control strategy that 
harnesses economies of scale. One approach is to use mi¬ 
crowave pulses from a single control source for multiple 
qubits [7], requiring frequency-matched qubits and high¬ 
speed routing of pulses to separate control lines. The lin¬ 
ear scaling of control equipment could then be reduced 
to a constant overhead for the most expensive resources. 

Using control equipment for multiple qubits has previ¬ 
ously been demonstrated for optical addressing in atomic 
systems, where qubits naturally have the same fre¬ 
quency [HHII]. Such frequency reuse also becomes pos¬ 
sible in circuit quantum electrodynamics (cQED) |12j in 
the context of fault-tolerant computation strategies [131- 
[H] which rely only on local interactions between qubits 
mediated by bus resonators. The natural isolation be¬ 
tween different lattice sites allows the use of repeating 
patterns of qubit frequencies with selectivity provided by 
spatial separation. A tileable unit cell with a handful of 
qubit frequencies m could therefore provide a promis¬ 
ing route towards scalability. Crucially, this also solves 
the frequency-crowding problem which arises when try¬ 
ing to fit many distinct-frequency qubits within the finite 
useful bandwith of the circuit-based devices, particularly 
for designs based on weakly anharmonic qubits where 
higher levels must also be avoided HadS]. While no 
qubit experiments have yet shown the viability of this 


approach, Hornibrook et al. have recently demonstrated 
a cryogenic switching matrix for pulse distribution op¬ 
erating at 20 mK, triggered by a field-programmable 
gate array at 4 K |7]. Cryogenic control equipment may 
shorten feedback latency and reduce wiring complexity 
across temperature stages, but the isolation and opera¬ 
tional frequency range achieved are currently insufficient 
for typical cQED experiments. 

In this Letter, we demonstrate frequency reuse in an 
extensible solid-state multiqubit architecture. Specif¬ 
ically, we show independent simultaneous control of 
two same-frequency qubits with a home-built room- 
temperature vector switch matrix (VSM). The VSM al¬ 
lows tailoring of control pulses to individual qubit prop¬ 
erties, and routing of the pulses to either one or both of 
the qubits using fast digital markers. We develop sev¬ 
eral different approaches to selective pulse broadcasting, 
including a simple scheme for implementing independent 
Clifford control on an arbitrary number of qubits with 
a constant overhead in time. The device for this ex¬ 
periment is designed to allow testing in a circuit with 
the correct connectivity of a relevant surface-code lat¬ 
tice da HO]. Using randomized benchmarking (RB), we 
show that all control schemes exceed the fidelity thresh¬ 
old for surface code and are dominated by qubit relax¬ 
ation. We also develop a method for measuring leakage 
to the second excited state directly within the context of 
RB [211 [22]. We characterize the limitations of our sys¬ 
tem and find no major obstacles to scaling up to larger 
implementations. 

Demonstrating frequency reuse in a context relevant 
for increasing system sizes requires two key elements: a 
method for distributing control pulses to multiple qubits 
using a single qubit-control source, and a multiqubit de¬ 
vice containing same-frequency qubits with relevant con¬ 
nectivity. Here, we focus on a particular implemen¬ 
tation of the surface code, where the connectivity be¬ 
tween nearest-neighbor qubits is achieved via bus res¬ 
onators m- Figure illustrates a conceptual design 
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FIG. 1. (color online). Schematic of independent control of 
same-frequency qubits using a vector switch matrix (VSM). 
The device (gold-colored box on right) connects two trans- 
mons with matched frequency (Qt and Qb: 6.220 GHz) indi¬ 
rectly through coupling buses and a third non-matched trans- 
mon (Qm: 6.550 GHz). This provides the smallest relevant 
subunit of the four-frequency surface-code fabric illustrated 
above right. The VSM (blue box) allows independent, si¬ 
multaneous transmon control with tailored DRAG pulsing by 
combining and directing Gaussian and derivative-of-Gaussian 
input pulses to separate outputs T and B connected to ded¬ 
icated drive lines for each qubit. The link between inputs 
and outputs can be switched on nanosecond timescales using 
the digital marker inputs Mt and Mb (orange lines). DRAG 
pulses for the targeted qubit are independently tuned in both 
amplitude and phase for each input-output pair (top left). 


based on repeated tiling of a unit cell consisting of four 
qubits with unique frequencies that are coupled via bus 
resonators. 

Our device contains a small block of this design, 
consisting of two same-frequency transmon qubits (Qt 
and Qb), which are connected to a third qubit (Qm) 
via separate bus resonators (Fig. [^. Each qubit has 
a capacitively-coupled drive line for individual qubit 
control [23], a readout resonator coupled to a com¬ 
mon feedline for frequency-division multiplexing read¬ 
out [231 [^, and a flux-bias line for individual fre¬ 
quency tuning (Fig.[^. While Qt and Qb were designed 
to be identical, fabrication uncertainties resulted in a 
sweet-spot (maximum) frequency of Qt 57 MHz higher 
than that of Qb. With Qb and Qm kept at their re¬ 
spective sweet-spots (6.220 GHz and 6.550 GHz, respec¬ 
tively), Qt was then flux tuned to match Qb with an 
accuracy of 50 kHz, determined using Ramsey measure¬ 
ments. The coherence times at the operating frequency 
can be found in the Supplemental Material [2^1 . Because 
of the transmon’s weak anharmonicity [27j . high-fidelity 
fast single-qubit control is achieved using the method 
of derivative-removal-via-adiabatic-gate (DRAG) puls¬ 
ing, where the in-phase Gaussian pulse is combined with 
an in-quadrature derivative-of-Gaussian pulse [28l l29] . 
For each qubit, this requires independent amplitude con¬ 
trol of the two constituent quadrature pulses. 


The VSM was designed to accept multiple input pulses 
and selectively fan them out to multiple qubits with in¬ 
dividual pulse tuning for each qubit (Fig. [^. Our home- 
built room-temperature 4x2 (four input, two output) 
VSM allows independent control of amplitude and phase 
for each of its input-output combinations. Fast marker- 
controlled digital switches enable routing of pulses to 
the qubits at nanosecond timescale, with approximately 
50 dB isolation in the frequency range from 4 to 8 GHz 
(see |26| for VSM specifications). By directing the two 
consistuent pulses of DRAG control through separate in¬ 
puts of the VSM, this allows independent, in-situ DRAG 
tuning for both same-frequency qubits using four AWG 
channels [26] . 

The first critical test of our control architecture is to 
assess the VSM’s ability to implement high-precision con¬ 
trol of one qubit while leaving the other qubit idle. To 
do this, we use the standard technique of single-qubit 
RB based on Glifford gates [30U32] . which allows us to 
characterize control performance independently of state 
preparation and measurement errors. After initializing 
all qubits in the ground state by relaxation, we use the 
VSM to selectively apply random sequences of Cliffords 
gates to only one of the same-frequency qubits, with 
lengths m ranging from 1 to 800 gates and the results 
averaged for 50 different random sequences. We decom¬ 
pose each gate into the standard minimal sequence of tt 
and ±7r/2 pulses around the x and y axes (16 ns pulses 
separated by a 4 ns buffer; tp = 20 ns) [22], requiring 
on average {Np) = 1.875 pulses per Clifford. This is in 
contrast to atomic pulses, where the 24 single-qubit Clif¬ 
fords can each be implemented with a single pulse [33j . 
After applying a final Clifford that inverts the cumula¬ 
tive effect of all m previous Cliffords, the driven qubit 
is ideally returned to the ground state, but as a result 
of imperfections such as gate errors and decoherence, 
the final ground-state population decays as a function 
of TO. The decay rate can be related to the average fi¬ 
delity per Clifford Fq [SO] [31] . In a strictly two-level sys¬ 
tem, the measured ground- and excited-state populations 
averaged over many sequences ((Pq) and (Pi)) both con¬ 
verge to 0.5 for large to. For weakly anharmonic trans¬ 
mon qubits, leakage to the second-excited state can be an 
important additional source of gate error, which can lead 
to a shift of the asymptotic value away from 0.5. We ad¬ 
dress this issue by performing the RB protocol both with 
and without an additional final tt pulse |34] , which allows 
us to explicitly estimate the populations of the first three 
transmon states (see ]2S] for details). 

We implement the above characterization for one qubit 
at a time by switching off the marker for the undriven 
qubit at the VSM [Fig. [^a,b)], in each case measuring 
the effect on both qubits simultaneously via multiplexed 
readout. From the results in Fig. ^c,d), we calculate 
the average Clifford fidelities for the two individually- 
driven qubits to be 0.9982(2) (Qt) and 0.9986(2) (Qb), in 
both cases surpassing the best known surface-code fault- 
tolerance threshold for single-qubit gates of ~ 0.99 [351- 
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FIG. 2. (color online). Single-qubit control of same-frequency 
qubits using the VSM. (a,b) Schematics showing DRAG 
pulses routed exclusively to either Qt or Qb using the corre¬ 
sponding markers (always on or off). Details of the tune-up 
procedure for optimizing pulse amplitude and phase parame¬ 
ters for each qubit are provided in Ref. 1261 (c,d) Gharacter- 
ization of single-qubit control by randomized benchmarking 
(RB) of Glifford gates. Average populations of Qt and Qb in 
the ground, first- and second-excited states {{Pq}, (Pi) and 
(P 2 ), respectively) as a function of the number of Glifford 
gates applied. Gurves are best fits of single exponentials with 
offsets. Single-qubit Glifford fidelities for each qubit are ex¬ 
tracted from the decay of ground-state populations. For both 
qubits, the fidelity surpasses the surface-code fault-tolerance 
threshold and lies close to the theoretical value expected for 
Ti-only relaxation-limited performance. (e,f) Expanded plots 
of second-excited state leakage during RB. Gurves are best fits 
according to Eq. <§• (g,h) Gross-excitation of the undriven 
qubit resulting from control pulses applied to the driven qubit. 
Further measurements indicate that this effect can be at¬ 
tributed to microwave leakage. 


ET]. We compare these values with the expected average 
Clifford fidelities assuming only Ti decay [38] : 
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This shows our results are predominately limited by re¬ 
laxation effects, the difference in performance being con¬ 
sistent with the different Ti times for the two qubits. 
Additional measurements in the Supplemental Material 
furthermore show that there is no difference in perfor¬ 
mance when both qubits are driven simultaneously by 
the same pulse sequence (both markers on). From the 
measured leakage populations (P 2 ), we also extract esti¬ 
mated per Clifford leakage rates k of 4.1(2) x 10“® (Qt) 
and 1.3(4) x 10“® (Qb) by fitting the following simple 


model to the data (see for details): 

(P2N) - (1 - , (2) 

where T 2 _>.i is the second- to first-excited-state relaxation 
time. As these leakage rates are much smaller than the 
gate errors (1—Fc), it is reasonable to neglect them when 
estimating the Clifford fidelity. 

We next explore the effect of the single-qubit control 
pulses on the undriven qubit [Fig. [^g,h)], which should 
ideally remain in the ground state. We fabricated the 
qubits as close to each other as possible in order to study 
the worst-case scenario for such cross-excitation effects. 
While Qt remains largely unaffected when driving Qb, 
a substantial deviation from the ground state is mea¬ 
sured in Qb when driving Qt. There are several possi¬ 
ble mechanisms for cross-excitation effects in our system: 
cross-coupling (higher-order quantum coupling mediated 
by the bus resonators and Qm) and cross-driving (spu¬ 
rious driving of the idle channel by microwave leakage 
resulting from imperfect isolation either in the VSM or 
on chip). From single-excitation swap experiments (see 
Ref. EH])) we observe a residual exchange interaction (T^j 
between Qt and Qb with strength J/27r < 36 ± 1 kHz. 
This symmetric swapping of excitation is therefore un¬ 
likely to explain the strong asymmetry in the amount of 
cross-excitation measured for the different qubits. Fur¬ 
thermore, in RB experiments, where the state of the 
driven qubit is moved randomly around the Bloch sphere 
during each run of the experiment, we expect the effect 
of cross-coupling to be dramatically reduced, because the 
period tt/J is far longer than the average Clifford gate 
time. Direct, independent measurements of microwave 
isolation in both the VSM and on chip indicate that the 
on-chip cross-driving significantly dominates [26] . This 
was also confirmed in-situ by performing the same RB 
measurements with the undriven qubit physically dis¬ 
connected from the VSM. This showed no significant 
reduction in cross-excitation. Furthermore, numerical 
simulations show that the observed effects are consistent 
with cross-driving alone at levels similar to the measured 
values [55]. We note, however, that while the plots in 
Fig.llg ,h) are useful diagnostics for identifying the pres¬ 
ence of the cross-driving effect, they should not be inter¬ 
preted in the same way as the RB plots for the driven 
qubit. A more comparable way to study this effect is to 
use interleaved RB [39| , where pulses on the driven qubit 
(ideally identity operations on the undriven qubit) are 
interleaved with a random sequence of Cliffords applied 
to the undriven qubit |26| . From this, we estimate the 
average idling fidelity for Qb to be 0.9986(5), which is 
consistent with the error due to Ti decay during idling. 
This confirms that cross-excitation effects do not domi¬ 
nate the error per Clifford, as characterized by RB. 

The defining test of extensibility in our control archi¬ 
tecture is to demonstrate the simultaneous, independent, 
single-qubit control over same-frequency qubits that is 
enabled by selective broadcasting using the VSM. We ex- 
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FIG. 3. (color online). Selective broadcasting schemes for 
simultaneous single-qubit control of multiple qubits. Main 
figure: Example of a single Clifford round for n=2 qubits 
targeting C 2 {Cis) in Qt (Qb) {Cj defined in Ref. l26ll . In 
the sequential scheme, the pulses implementing C 2 are di¬ 
rected to Qt, after which the pulses implementing Cia are 
directed to Qb . In the compiled scheme, the two Cliffords are 
realized concurrently using a pre-determined pulse sequence, 
with appropriate markers, which minimizes the total number 
of pulses, A^p (see [26] for the compilation algorithm). Finally, 
in the 5-primitives scheme, a fixed sequence of five pulses is re¬ 
peated in each round (Afp=5). The targeted Cliffords are then 
applied simultaneously by selecting the appropriate subset of 
pulses for each qubit (see [^ for the 5-primitives marker ta¬ 
ble). Top-right: Scaling of the average pulses per multiqubit 
combination of Cliffords, {Np), versus qubit number n. The 
constant scaling achieved by the 5-primitives scheme provides 
a dramatic improvement over the linear scaling of the sequen¬ 
tial scheme. While {Np) is always lowest for the compiled 
scheme, pre-compiling the optimal pulse and marker combi¬ 
nations is impractical for n > 5, and the improvement over 
the simpler 5-primitives scheme is negligible by n ~ 10. 
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FIG. 4. (color online). Characterization of sequential, com¬ 
piled and 5-primitives selective-broadcasting schemes by RB. 
(a,b) Evolution of the average transmon populations for Qt 
(a) and Qb (b) as a function of the number of Clifford rounds. 
Curves are the best fits of single exponentials with offsets. 
(c,d) Expanded plots of second-excited-state leakage during 
RB. (e) Average single-qubit Clifford-gate fidelities for Qt 
and Qb in each scheme, extracted from the decay of the cor¬ 
responding ground-state population. All fidelities surpass the 
surface-code fault-tolerance threshold and closely track those 
expected for Ti-relaxation-limited performance. Compiled se¬ 
lective broadcasting performs best, as expected, having the 
lowest total number of pulses (see Fig. [^. 


plore three paradigmatic schemes for implementing selec¬ 
tive broadcasting of Cliffords on an arbitrary number of 
qubits n (Fig. [^. In the most straightforward selective- 
broadcasting scheme, the individual qubits are driven se¬ 
quentially, with each pulse being directed to one qubit 
at a time. This results in a linear scaling of the average 
number of pulses per Clifford round ((IVp) = 1.875 x n). 
By contrast, the second paradigm takes best advantage of 
the VSM’s capability to broadcast simultaneously to mul¬ 
tiple qubits by compiling the constituent Clifford pulses 
to minimize Np for each Clifford combination in the se¬ 
quence. The compilation is performed by searching all 
possible combinations of single-qubit Clifford decompo¬ 
sitions and finding the one that minimizes Np (see [25] for 
further information). In Fig.j^ we show exact values and 
estimates for {Np) for the compiled scheme up to n = 10. 
Unfortunately, the compilation run-time increases expo¬ 
nentially with the number of qubits. This motivates our 
final broadcasting paradigm, where all Clifford gates can 
be implemented using the same fixed, ordered sequence of 
five pulse primitives (Fig. [^. Independent Cliffords can 
be applied to all qubits, irrespective of n, by selectively 


directing the appropriate subset of pulses to each qubit, 
achieving a constant overhead in time for single-qubit 
Clifford control. While the compiled scheme by defini¬ 
tion always provides the minimum sequence length, our 
estimates of {Np) for the compiled scheme suggest that 
it asymptotes to the same value of 5 achieved by the 
simple, prescriptive 5-primitives scheme. We note that 
the specific choice of pulse primitives is not unique, but 
at least five primitives are required (four pulses allow a 
maximum of 16 unique gate decompositions, compared 
with the 24 single-qubit Cliffords). 

To demonstrate the full functionality of this control ar¬ 
chitecture, we implement all three selective-broadcasting 
schemes and measure their performance using parallel 
single-qubit RB with independent Clifford sequences for 
each qubit. Figure |^ shows that the compiled scheme 
performs best, followed by the sequential and then 5- 
primitives schemes, consistent with the average number 
of pulses required for each (Fig. |^. In all cases, the av¬ 
erage fidelity per Clifford still surpasses the surface-code 
fault-tolerance threshold, and the average error is again 
dominated by relaxation. The results are completely con- 
































5 


sistent with the values obtained in the test for isolated 
single-qubit control, indicating no substantial decrease in 
gate performance using selective broadcasting schemes. 

Our VSM allows efficient reuse of qubit control equip¬ 
ment on same-frequency qubits. It enables high-precision 
single-qubit control of multiple qubits with a perfor¬ 
mance that is mainly limited by relaxation. We have 
demonstrated three selective broadcasting schemes, all 
of which achieve a performance that surpasses the fault- 
tolerance threshold for the surface code. In particular, 
the 5-primitives scheme implements arbitrary Clifford 
control with a fixed five-pulse sequence, where the target 
Clifford is selected by routing a subset of the pulses using 
digital markers. By adding a sixth, non-Clifford gate to 
the five pulse primitives, this can be extended to achieve 
universal single-qubit control. Combining the connectiv¬ 
ity of our device, the VSM-based control, and the hxed 
pulse overhead of the 5-primitives broadcasting strategy, 
our experiment realizes the simplest element of an exten¬ 


sible qubit control architecture. While we do not yet see 
explicit savings in control hardware for two qubits, this 
design can be expanded to more same-frequency qubits 
without any further increase in microwave sources or ar¬ 
bitrary waveform generators. This experiment suggests 
that surface-code tiling with frequency reuse is a viable 
path towards large-scale quantum processors. 
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SUPPLEMENTAL MATERIAL 

This supplement provides experimental details and ad¬ 
ditional data supporting the claims in the main text. 
The device and experimental setup, including images 
of the device and a full wiring diagram, are described 
in Section |T] The microwave performance of the vec¬ 
tor switch matrix (VSM) and its use for independent 
control of two same-frequency qubits are demonstrated 
in Sec. |II A| The techniques employed for tuning qubit 
pulses are discussed in Sec. Ell Section |TV] contains the 
results of randomized benchmarking (RB) experiments 
in the two-qubit global broadcasting context. Our tech¬ 
nique for assessing the effects of leakage is detailed in 
Sec. |Vj Cross-coupling and cross-driving effects are char¬ 
acterized in Sec. lYll along with numerical simulations of 
the effects of cross-excitations on RB. In Section [Vl D[ we 
introduce a method for generating selective-broadcasting 
pulse sequences that are robust to cross-excitation effects. 
The decompositions of the 24 single-qubit Cliffords into 
a minimal set of pulses and the 5-primitive pulses are 
given in Sec. |VII| Finally, the algorithm used to compile 
optimal pulse sequences for implementing independent 
single-qubit Clifford gates on multiple qubits is explained 
in Sec. IVIIIl 


I. QUANTUM CHIP AND EXPERIMENTAL 
SETUP 

A. Chip design and fabrication 

Our quantum chip consists of three transmons 
(top: Qt, middle (ancilla): Qm, and bottom: Qb) with 
dedicated voltage drive lines (Dt, Dm and Db, respec¬ 
tively), flux-bias lines, and readout resonators. All read¬ 
out resonators are capacitively coupled to one common 
feedline which crosses various on-chip components us¬ 
ing airbridge crossovers. Qubits Qt and Qm are cou¬ 
pled by one bus resonator, and Qm and Qb by another 
(fundamental frequencies 4.9 and 5.0 GHz, respectively). 
All resonators are open-ended on the coupling side, and 
short-circuited at the other. 

The chip fabrication method is similar to that in Ref.[l] 
but with some important differences which we now ex¬ 
plain. Rather than sapphire, we use a high-resistivity in¬ 
trinsic silicon substrate prepared by HF dip and HMDS 
surface passivation before sputtering a 300-nm-thick film 
of NbTiN, as introduced in Ref. [51 This change aims 
to improve the substrate-metal interface and thereby in¬ 
trinsic quality factors for both resonators and qubits. 
After sputtering, the patterns are etched into the su¬ 
perconducting layer using reactive-ion etching with an 
SF 6 /O 2 plasma. In contrast with the A1 transmon ca¬ 
pacitor plates commonly used, in this experiment we 
make them also from NbTiN, with an aim to improve 
the substrate-metal interface and avoid large AlO^, sur¬ 
faces which may house unwanted two-level systems. Only 



feedline Qj, Qb and drive lines 

readout resonators Qm and drive line 

bus resonators flux-bias lines 


FIG. SI. Scanning electron micrograph (with added false 
colour) of a device twin to the one used in the experiment and 
fabricated in the same batch. Notably, the feedline crosses the 
three readout resonators. 


the Josephson junctions are made by the standard tech¬ 
nique of A1-A102,-A1 double-angle evaporation. A further 
HF dip just prior to evaporation also helps to contact 
the junctions directly to the NbTiN capacitor plates. In 
Ref. [U air bridges were already used to cross the feed¬ 
line over flux-bias lines on chip. Here, we extend this 
technique to allow the feedline to cross three readout res¬ 
onators. 

A key requirement for this experiment was the ability 
to match qubit frequencies without sacrihcing coherence. 
Flux-bias lines allow easy compensation for mismatch, 
but at the cost of reduced coherence in the qubit detuned 
from its maximal frequency (coherence sweet spot [3]). 
We aimed for identical maximum frequencies of Qt and 
Qb, as determined by capacitor and junction geometries. 
The capacitors were easily matched in fabrication. We 
then selected the chip with the closest matching room- 
temperature resistance values for the relevant qubit junc¬ 
tions. 


B. Experimental Setup 

The chip is anchored to a copper cold-finger connected 
to the mixing chamber of a Leiden Cryogenics CS8I 
^He/'^He dilution refrigerator with 7 mK base temper¬ 
ature. A copper can seals the sample space, with an 
inner surface that is coated with a mixture of Stycast 
2850 and silicon carbide granules (15 to 1000 nm diam¬ 
eter) used for infrared absorption [4]. The copper can is 
in turn magnetically shielded by an aluminum enclosure 
and two outer Cryophy enclosures (1 mm thick) [5]. 

A complete wiring schematic showing all cryogenic 












and room-temperature components is shown in Fig. S2 


The four analog channels of the Tektronix AWG5014C 
create the in-phase and in-quadrature pulses for Qx 
and Qb by single-sideband modulation of a common 
carrier. Because single-sideband modulation requires 
two AWG channels to modulate an IQ mixer, indepen¬ 
dent derivative-removal-via-adiabatic-gate (DRAG) tun¬ 
ing with the VSM therefore requires four AWG channels, 
irrespective of the number of qubits. The VSM can be 
scaled up to many output channels, and direct hardware 
savings can be realized as soon as three or more same- 
frequency qubits are driven by a single set of AWG in¬ 
puts. These pulses are input at ports 1 and 2 of the VSM. 
The VSM combines these pulses with individually tuned 
insertion loss and phase to each of two outputs (labelled 
T and B). Input-output combinations can be switched on 
nanosecond timescales using the gate inputs of the VSM, 
provided by digital markers of the AWG5014G. A sec¬ 
ond AWG5014C with the appropriate carrier frequency 
is used to excite transmons Qx and Qb to the second- 
excited state (Sec. |V]), to pulse measurement tones, and 
to trigger the AlazarTech ATS9870 acquisition card. 



FIG. S2. Detailed schematic of the experimental setup and 
optical image of the chip (installed in a Leiden Cryogenics 
CS81 dilution refrigerator). Note that while our VSM has two 
pairs of analog inputs [(1,2) and (3,4)] and four gate inputs 
(activating links from each pair to each output), the pair (3,4) 
and its associated gate inputs are not used throughout this 
experiment. 
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C. Device frequencies 


Qubit 

/,nax (GHz) 

/bias (GHz) 

/res (GHz) 

Qt 

6.277 

6.220 

6.700 

Qm 

6.551 

6.551 

6.733 

Qb 

6.220 

6.220 

6.800 


TABLE SI. Table of sweet-spot frequencies /max and bias- 
point frequencies /bias of the three qubits, as well as the funda¬ 
mental frequencies /res of their dedicated readout resonators 
at the bias point. The qubits are tuned into the bias point by 
a combination of spectroscopy and standard Ramsey experi¬ 
ments. 


D. Qubit coherence times 



Frequency (GFIz) 



Qt 

Qm 

Qb 

1.00 

^ =8.7 fis 

^ =7.4 fis 

< ' =11.1^5 ■ 

0.75 




0.50 




0.25 




0.00 





0 5 10 15 20 

6 5 10 15 20 

0 5 10 15 20 

1.00 

i T 2 —4.7 fis 

> ' ' T2*=9.7'irs 

S n ^2* =12.5 fis' 

0.75 




cC 0.50 




0.25 




0.00 





0 2 4 6 8 10 

6 2 4 6 8 10 

0 2 4 6 8 10 

1.00 

T 2 ^8.9 fis ■ 

' ' ' =14.8 IIS 

' ' /f bi7,b iis 

0.75 

f \ - 

A A A/\. 

A Aa/v 

0.25 


/ v V 

:/ V V v\; 

0.00 

1 .- 

/. , , , . 

/ .- 


0 5 10 15 20 25 0 5 10 15 20 25 0 5 10 15 20 25 

Wait time {iis) 


FIG. S3. Measurements of relaxation (Ti, top), Ramsey de¬ 
phasing {T 2 , middle) and echo dephasing (T®, bottom) times 
for the three qubits at the bias point. When measuring Qt 
or Qb, the other qubit is detuned by —50 MHz to suppress 
cross-coupling effects. Ramsey fringes for Qt (middle, left 
panel) fit better to a Gaussian (shown) than an exponential 
decay, reflecting the susceptibility of Qt to low-frequency flux 
noise away from its sweet-spot. Pi denotes excited-state pop¬ 
ulation. 


FIG. S4. Insertion loss between inputs and outputs of the 
VSM for four static combinations of the gate inputs. The in¬ 
sertion loss is measured relative to the level with both gate 
inputs activated (outoub). The black curves indicate the 
noise background in our scalar network analyzer measure¬ 
ment. The dashed vertical lines indicate the common fre¬ 
quency (6.220 GHz) of Qt and Qb at the bias point. 


B. Individual qubit tune-up 


The VSM enables independent control of the on/off 
state, insertion loss and phase for every input-output 
combination. We exploit this feature to perform DRAG- 
compensated pulses on Qt and Qb, that are individually 
tailored for each qubit. Different types of gate errors, 
such as non-ideal in-phase and in-quadrature amplitudes, 
can be distinguished using an AllXY sequence [S], con¬ 
sisting of 21 combinations of two puls es drawn fro m the 
set (Table] 


S2|. Figure 


S5 


shows 

AllXY sequence results for Qt and Qb as the amplitude 
of each quadrature on Qt is varied independently. While 
the AllXY signature of Qt reveals changing levels of am¬ 
plitude and phase errors, there is no noticeable change in 
the AllXY signature of Qb. This demonstrates the use 
of the VSM for individual tune-up of pulses for same- 
frequency qubits. 


II. VECTOR SWITCH MATRIX 


A. Measured isolation 


To characterize the isolation of the VSM in the range 
4 to 8 GHz, we have measured the insertion loss between 
all input (1 and 2) and output ports (T and B) with static 
settings at the two gate inputs (Fig. S4). Ideally, each 
gate activates (on state) and deactivates (off state) the 
link of both inputs 1 and 2 to one output, independent of 
the other gate. As shown in Fig. |S4[ the typical relative 
isolation with the relevant gate in the off state is ~ 50 dB. 
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(a) Qx (b) Qb 




AlIXY ID 


FIG. S5. AllXY results for Qt (left) and Qb (right), as the in- 
phase amplitude (top) or in-quadrature amplitude (bottom) 
of Qt pulses is varied. During each measurement, both qubits 
are driven simultaneously, but read out individually. The 
AllXY results of Qb all overlap independent of the setting on 
Qa, showing independent control of pulse parameters for each 
qubit. 


ID 

Ideal Pi 

Pulses 

First Second 

ID 

Ideal Pi 

Pulses 

First Second 

1 

0 

I 

I 

12 

0.5 

Y„ 

U/2 

2 

0 


x^ 

13 

0.5 

x^ 

X^/2 

3 

0 

Y-K 

Y^ 

14 

0.5 

X„/2 

x„ 

4 

0 


Y^ 

15 

0.5 

X-r, 

X„/2 

5 

0 

Y-K 

x^ 

16 

0.5 

U/2 

Y^ 

6 

0.5 

I 

Xti-/2 

17 

0.5 

Y„ 

to 

7 

0.5 

I 

U/2 

IS 

1 

I 


8 

0.5 

X^/2 

U/2 

19 

1 

I 

Y^ 

9 

0.5 

U/2 

X„/2 

20 

1 

X^/2 

X„/2 

10 

0.5 

X^/2 

Y„ 

21 

1 

U/2 

U/2 

11 

0.5 

U/2 

x^ 




TABLE S2. The 21 two-pulse combinations comprising the 
AllXY pulse sequence [5]. 



FIG. S6. Fine calibration of pulse amplitude by initial 7r/2 
pulse, followed by 2N repeated tt pulses. The initial slope de¬ 
termines if the qubit is over- (positive slope) or under-driven 
(negative slope). 


five or negative initial slope on Pi versus iV, respectively 
(Fig. S6). We choose the in-quadrature amplitude that 
minimizes the absolute slope. 


B. DRAG-parameter calibration 

To minimize phase errors resulting from the presence 
of the second- and higher-excited states, we optimize 
the scaling of the in-quadrature pulse. As in conven¬ 
tional DRAG [H [7], we choose as the envelope of the 
in-quadrature pulse the derivative of the Gaussian enve¬ 
lope on the in-phase pulse. The DRAG scaling param¬ 
eter is calibrated using the method detailed in Ref. [5]. 
Specihcally, we measure the difference in excited-state 
population produced by the and XT^Y^f 2 pulse 

combinations (AllXY ID 10 and 11). Ideally, for both, 
the final qubit state would lie on the equator. However, 
any phase error shifts the final excited-state population in 
opposite directions in these cases. We choose the DRAG 
scaling parameter minimizing this shift. 


IV. GLOBAL BROADCASTING 


III. PULSE-CALIBRATION ROUTINES 

We tune up qubit pulses by alternating the calibra¬ 
tion of in-phase and in-quadrature pulse amplitudes until 
a simultaneous optimum is found. The two calibration 
routines are discussed below. 


A. Accurate in-phase pulse amplitude calibration 

The in-phase quadrature amplitude is calibrated by 
first applying a 7r/2 pulse to the qubit, followed by a 
train of tt pulses. The pulse sequence is (A.^.) Xt ^/2 |g), 

where \g) is the qubit ground state and N S [0,49]. In 
the absence of gate errors and decoherence, the driven 
qubit would end on the equator of its Bloch sphere for 
all N. However, over- or under-driving produces a posi- 


Aside from single-qubit control and selective broad¬ 
casting, the VSM also allows global broadcasting of 
pulses to all qubits simultaneously by keeping the mark¬ 
ers for both qubits on (Fig. [S^ . While this does pro¬ 
vide simultaneous control of Qt and Qb, marker control 
is needed to achieve independent control. Using RB, we 
measure the performance of both qubits when broadcast¬ 
ing pulses to both qubits, and compare the results with 
those obtained from single-qubit control (Fig. S7). The 
global broadcasting RB measurements were alternated 
with the single-qubit RB measurements, and aside from 
marker settings all other settings were identical. Gom- 
parison of the results in Fig. ^ show that the qubit gate 
performance does not depend on whether a single qubit is 
controlled, or both are controlled simultaneously through 
global broadcasting. 
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Number of Clifford gates m 


FIG. S7. Global broadcasting of DRAG pulses to same- 
frequency qubits, (a) Illustration of global broadcasting. Two 
simultaneous pulses, one with Gaussian envelope at input 1, 
and another with derivative-of-Gaussian envelope at input 2, 
are simultaneously directed to Qt and Qb (both markers al¬ 
ways on). The insertion loss and phase shift of each pulse 
is separately optimized for each output to produce precision 
DRAG pulses for each qubit. (b,c) Gomparison of single-qubit 
driving versus driving both qubits (broadcasting) by RB of 
Clifford gates composed from tt/2 and tt pulses [8]. Average 
population of Qt and Qb in the ground, first- and second- 
excited state {{Po}, (Pi) and {P 2 }, resp.) as a function of 
the number of Clifford gates applied. Curves are the best fits 
of single exponentials with offsets to the populations. The 
single-qubit Clifford-gate fidelity for each qubit is extracted 
from the decay of the corresponding ground-state population 
when using global broadcasting. 


V. LEAKAGE TO SECOND EXCITED STATE 


assume that leakage is irreversible. We therefore model 
leakage using the following difference equation for {P 2 )' 

{P2[m+1]) - {P2[m]) - tp{Np) K-{P2[m]), (S2) 

J-2^1 

where T 2 _j.i is the second- to first-excited-state relax¬ 
ation time. Assuming no initial population in the second- 
excited state, the solution is Eq. (2), which shows good 
agreement with measured data. We extract k by fitting 
Eq. (2) to (P 2 ) data. T 2 ^i is obtained from the best-fit 
decay constant (not directly measured) and k from the 
best-fit prefactor. 


VI. CROSS-COUPLING AND CROSS-DRIVING 
EFFECTS 

There are several sources of spurious cross-qubit inter¬ 
actions [see Fig. 2(g,h) of the main text] which play a 
role in our experiments. We divide these into two main 
classes: cross-coupling, where the qubits themselves are 
coupled via a direct or indirect quantum interaction, and 
cross-driving, where input microwaves directly drive the 
untargeted qubit with a reduced amplitude as a result 
of imperfect isolation either on or off chip. These cross¬ 
excitation effects depend strongly on specific chip design, 
in our case chosen according to our primary aim to study 
the potential of frequency reuse in a circuit fully compati¬ 
ble with a larger surface-code lattice. This governed both 
how qubits and resonators were connected, as well as the 
selection and arrangement of frequencies. In addition, 
in order to fully explore the limitations of such tech¬ 
niques, the qubits were positioned on the chip as close 
to each other as possible to provide a worst-case scenario 
for cross-excitation effects. 


Leakage is fundamentally different from unitary qubit 
errors. To quantify leakage, we monitor the populations 
Pi of the three lowest energy states {i G {0,1, 2}) during 
RB and calculate the average values {Pi) over all seeds. 
To do this, we calibrate the average signal levels T) for the 
transmons in level i, and perform each RB measurement 
twice, the second time with an added final tt pulse on 
the 0-1 transition. This final tt pulse swaps Pq and Pi, 
leaving P 2 unaffected. Under the assumption that higher 
levels are unpopulated {Pq + Pi + P 2 = 1), 


V0-U2 VI-V 2 

>o' 


'S-V 2 ' 

Vi — V 2 Vo — V 2 

Pi 


S' - V2_ 


(SI) 


where S {S') is the measured signal level without (with) 
final TT pulse. The populations are extracted by matrix 
inversion. 

Measuring (P 2 ) as a function of the number of Clifford 
gates allows us to estimate an average leakage per Clif¬ 
ford, K. Because the populations are ensemble averages 
over different random seeds, we assume that leakage of 
the average qubit-space populations to {P 2 ) is incoher¬ 
ent, and, provided {P 2 ) remains small {k small), we also 


A. Cross-coupling 


In our device, the same-frequency qubits Qt and Qb 
are connected through a linear chain of coupled quan¬ 
tum elements: the two bus resonators and intermedi¬ 
ate ancilla qubit Qm. They are also coupled capaci- 
tively through the ground plane. When the intermedi¬ 
ate elements are detuned in frequency away from the 
near-resonant qubits, such geometries typically give rise 
to higher-order exchange-type interactions between the 
qubits of the form J{a^a^ -I- o't'^b) [Sllin]- Consistent 
with this picture, we observe coherent swapping of exci¬ 
tations between Qt and Qb at a rate strongly dependent 
on the qubit detuning (Fig. S8). 


In order to characterise the cross-coupling between Qt 
and Qb, we measure the evolution of excited-state pop¬ 
ulations after a single excitation is injected at one of the 
qubits with a tt pulse. To place a tight upper bound on 
the interaction strength J, the qubit frequencies must be 
matched as closely as possible. We achieve an accuracy 
of around 50 kHz using Ramsey experiments, limited by 
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(a) Initially exciting Q,[, (b) Initially exciting Qg 



FIG. S8. Temporal evolution of a single excitation after ini¬ 
tially exciting Qt (a) or Qb (b). The oscillations of population 
in both qubits are out of phase and have a common frequency. 
The sum of both populations shows approximately exponen¬ 
tial decay, with time constant 10.2(2) /rs, intermediate be¬ 
tween the relaxation times of Qt and Qb. 




Normalized pulse amplitude 



o.oo l-. 

6.00 6.05 6.10 6.15 6.20 6.25 
Qj. frequency (GHz) 


(d) 


Qubit 

Driving strength 
Dt Db 

Qt 

1.0000 

0.0023 
(-53 dB) 

Qb 

0.0057 
(-45 dB) 

1.0000 


a combination of factors: the resolution of the flux tun¬ 
ing, the fitting resolution limit imposed by qubit T 2 de¬ 
phasing times, and also the frequency shifting induced 
by the qubit-qubit exchange interaction itself. As shown 
in Fig. S8 the total excitation number exhibits a typical 
exponential Ti decay, while the individual populations 
show a symmetric oscillation of the excitation between 
the two qubits. The measured common frequency of de¬ 
caying oscillations for both qubits sets an upper bound 
on the coupling strength of J/27r < 36 ± 1 kHz. 


B. Cross-driving 

When driving one qubit, the signal applied to its ded¬ 
icated voltage line residually drives the distant qubits 
on the chip. This results from microwave leakage in the 
VSM, near the sample, or from direct on-chip coupling 
of the drive line to the distant qubits. In Section p^I A[ we 
fully characterise the isolation in the VSM. This leads 
to leakage of around —57 dB and —54 dB on Qt and 
Qb, respectively, at the qubit operating frequency for 
the conditions used in the main experiments. To char¬ 
acterize the residual cross-driving at the device, we dis¬ 
connect the VSM and compare the amplitude required 
for pulses applied on one drive line (either Dt or Db) to 
perform a tt rotation on both Qt and Qb . We define the 
cross-driving strength of each drive line as the ratio Vc 
of the TT-pulse amplitude for the directly-driven to that 
for the cross-driven qubit. For this test, pulses are first 
amplified and then attenuated using a step attenuator to 
allow the large amplitude range required. Results shown 
in Fig. |S9| demonstrate on-chip cross-driving strengths 
below 1% (—53 dB and —45 dB on Qt and Qb, respec¬ 
tively), but still somewhat larger than cross-driving ef¬ 
fects resulting from finite isolation in the VSM. During 
the main experiments, both VSM and on-chip leakage are 
able to contribute to cross-driving effects. 

In the excitation-swap experiments of the previous sec¬ 
tion, it is clear the effects result from cross-coupling 
rather than cross-driving, because the exchange dynam- 


FIG. S9. Rabi oscillations of Qt (a) and Qb (b) induced by 
pulses on Dt at the bias point. The pulse amplitude is nor¬ 
malized to the TT-pulse amplitude for Qt. (c) Resonant cross¬ 
driving of Qb via Dt for a range of Qt frequencies. The weak 
dependence of the cross-driving strength shows that this ef¬ 
fect is most likely due to direct coupling of Dt to Qb. (d) 
Summary of cross-driving for Dt and Db at the bias point. 


ics occur while the drive pulses are off. Similarly, we can 
rule out these residual driving effects being related to 
cross-coupling, because the 16-ns Rabi pulses take only 
a fraction of the time required for a cross-coupling ex¬ 
change oscillation. As a further check that the effects 
observed are indeed due to cross-driving, we measure the 
cross-driving on Qb through Dt as Qt is tuned through 
resonance with Qb [Fig. [S^c)]. The cross-driving ratio 
is essentially constant, showing no significant dependence 
on the detuning between Qt and Qb. 


C. Cross-excitation in randomized benchmarking 


The isolated single-qubit control experiments in the 
main text [Fig. 2(g,h)] show that significant spurious ex¬ 
citations can build up in the idling qubit over the course 
of the long gate sequences tested in RB (particularly in 
the case of idling Qb while driving Qt). It may there¬ 
fore initially be somewhat surprising that virtually the 
same individual qubit-control performance is achieved 
in both selective-broadcasting (Fig. 4, main text) and 
global-broadcasting (Fig. S7| scenarios. As discussed in 


the main text, the observed cross-excitation is unlikely 
to result from cross-coupling, primarily because a sym¬ 
metrical quantum coupling should not result in strongly 
asymmetric effects on the different qubits. We now show 
the results are, however, consistent with the effects of 
cross-driving by numerically simulating RB with cross¬ 
driving under experimentally realistic conditions (using 
independently measured qubit and cross-driving param¬ 
eters). Simulations are performed using QuTiP [TTIIT^ . 
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Qt simulation Qg simulation 



Number of Clifford rounds in 

FIG. SIO. Simulation of cross-driving effects during RB. 
The results shown are averaged over ten runs, nsing the 
single-qubit minimal-set decomposition (red), and using 
the selective-broadcasting asymmetric and symmetric 5- 
primitives schemes (green and blue, respectively). Cross¬ 
driving effects are largely suppressed in the 5-primitives 
schemes by choosing the five pulse primitives such that con¬ 
stituent pulses largely cancel out. The symmetric 5-primitives 
scheme further reduces cross-driving effects by alternating be¬ 
tween the five pnlse primitives and the inverse pnlses. 


We model our system as two uncoupled qubits, Qt 
and Qb, subject to Ti relaxation (with corresponding re¬ 
laxation times) and cross-driving. We approximate the 
system dynamics using instantaneous unitary pulse oper¬ 
ators from the standard Pauli set W-, h±^/ 2 } 

with 20 ns delays of Ti-only qubit relaxation between 
pulses implemented using a master equation. When ap¬ 
plying a pulse to one qubit, cross-driving of the other 
qubit is implemented by applying a pulse with the same 
rotation axis, but with the original rotation angle mul¬ 
tiplied by the relevant cross-driving ratio. We note that 
also trying to model the effects of qubit dephasing using 
a simple master equation does not produce RB data con¬ 
sistent with the experimental observations (e.g., Figs 2 
and 4 of the main text). This reflects the non-uniform 
phase noise spectrum which affects the transmon qubit. 
The long RB pulse sequences consisting of tt and 7r/2 
pulses around the X and Y axes seem to provide some 
form of dynamical decoupling which makes the RB mea¬ 
surements robust to qubit dephasing. As will be seen, 
the experimental results can be well modelled using only 
Ti-type noise processes. 

Randomized benchmarking is implemented by gen¬ 
erating independent Clifford sequences for each qubit. 
We decompose Clifford gates using either the minimal 
set decomposition or one of the selective-broadcasting 
schemes. Figure |S10| shows simulated results of cross¬ 
driving for the isolated single-qubit control scenario re¬ 
ported in Fig. 2 of the main text. In this section, we are 
only concerned with the red curves, which correspond to 
implementing single-qubit RB with the standard set of 
pulse decompositions [5] . These simulations can be com¬ 
pared directly with the curves in Fig. 2(g,h). While the 
maximum excitation population observed in the simula¬ 
tions is larger than the value observed in the experiments, 
the simulations for both qubits show the same qualita¬ 
tive behaviour as the measured data. The quantitative 



Number of Clifford rounds in 


FIG. Sll. Simulation of sequential (interleaved) RB for sev¬ 
eral cross-driving ratios. For each cross-driving ratio, 50 
simulation runs were performed, using sequential RB up to 
m = 800 Glifford rounds. Under the cross-driving levels per¬ 
taining in our experiments (0.0037 and 0.0076 for Qx and Qb, 
respectively), the error per Glifford for the idling operation is 
dominated by the effects of Ti relaxation as calculated from 
Eq. (1). 


difference may be explained by the fact that the direct 
measurements of cross-driving were made at a different 
time from the main measurement run and we observed 
some small fluctuations in cross-driving levels over time. 

As discussed in the main text, while the plots of cross¬ 
excitation during RB are useful diagnostics of the pres¬ 
ence of a spurious cross-driving effect, they may give a 
misleading impression when presented in parallel with 
RB results. Although the decay curves look superfi¬ 
cially similar, they should not be interpreted in the same 
way. By contrast, the technique of interleaved RB (IRB), 
which was introduced to enable rigorous quantification of 
the performance of individual gates, allows us to calcu¬ 
late a meaningful error per Clifford for the idling opera¬ 
tion [13]. In IRB, the usual random sequence of Cliffords 
is alternated with identical repetitions of an individual 
gate. By comparing the interleaved decay rate with the 
decay rate for a standard RB measurement, it is possi¬ 
ble to calculate a robust error per gate for the individ¬ 
ual gate in question. In this context, the target gate is 
the nominal identity operation on one qubit which re¬ 
sults from a random Clifford being applied to the other 
qubit. The IRB pulse sequence is therefore identical 
to the sequence implemented in the sequential selective- 
broadcasting scheme (see Fig. 4 of the main text). In the 
main text, we use the formulas in Ref. m to calculate 
the idling error per Clifford, but for these simulations, 
the performance of sequential selective-broadcasting al¬ 
ready provides a simple way to assess the performance 
of cross-excitation during idling. Figure [STT] shows that 
idling performance as quantified by RB is limited mainly 
by Ti relaxation. Finally, when identical gate sequences 
are being applied to both qubits, cross-driving will result 
in a small amount of over-driving on each qubit (over¬ 
driving ratio To), which would also look like an error in 
pulse rotation angle. Figure [ST^ shows that the Clifford 
error is insensitive to both cross-driving and over-driving 
to first order. 





















14 



Rotation error 


Qt 


Qb 



FIG. S12. Simulations of Clifford fidelity Fc as a function of 
the cross-driving ratio (ret red) and the relative over-driving 
rotation error (ro —1: green). For each error type, 50 sim¬ 
ulation runs were performed, using sequential RB on both 
qubits simultaneously up to m = 800 Clifford rounds, after 
which Fc was extracted from an exponential fit to averaged 
data. The Ti limit [Eq. (1)] is given by a horizontal dashed 
line. The data shows that Fc is first-order insensitive to both 
cross-driving and to over-rotations. 


FIG. S13. Measured cross-driving results when performing 
selective broadcasting. Even though the selective broadcast¬ 
ing schemes are meant for multiple qubits, the markers of 
the measured qubit are turned off to measure cross-driving 
effects. The cross-driving effects of Qb are stronger than 
for Qt, in agreement with single-qubit results. With pulse 
decompositions whose cumulative effect largely cancels out, 
cross-driving is strongly reduced in the 5-primitives method. 


D. Making pulse sequences robust to cross-driving 

We have already shown that cross-excitation does not 
have a dominant effect on single-qubit control in both 
global and selective broadcasting. We show here that 
any residual effect can be largely eliminated also while 
a qubit is idling by choosing robust pulse sequences for 
decomposing the Clifford gates. 

If a qubit is idle, every pulse that is applied to the 
driven qubit rotates the idle qubit by an amount de¬ 
pending on the cross-driving ratio. The random applica¬ 
tion of successive pulses to the driven qubit can there¬ 
fore be viewed as a random walk for the idle qubit. 
As we will discuss in more detail in Sec. lYml there 
are many ways to compose a given Clifford gate from 
a small set of standard rotations. By choosing the con¬ 
stituent pulses in such a way that their combined ap¬ 
plication largely cancels out, cross-driving effects can be 
greatly reduced. In the standard set of Clifford decom¬ 
positions [8], the decompositions involve a majority of 
pulses rotating in the positive direction, biasing the ran¬ 
dom walk and producing a pronounced net cross-driving 
effect. This effect can be countered by choosing de¬ 
compositions which minimize the bias. We have imple¬ 
mented this in the 5-primitives scheme, by choosing the 
first three pulses, {Al,r/ 2 j^ 7 r/ 2 :-^ 7 r/ 2 }j to be positive ro¬ 
tations, and the last two, to move in the 

negative direction. Even though the pulse subset that is 
applied depends on the Clifford chosen, the pulses still 
largely cancel out after applying many Cliffords. Fur¬ 
thermore, as the single-qubit Clifford operations form a 
group, the inverse of all Cliffords also form the Clifford 
group. The complete inverse of the five pulse primitives, 
W, X_^/ 2 ,Y-k/ 2 ,X_^/ 2 }, can therefore also gener¬ 
ate each of the 24 Cliffords using an appropriate sub¬ 
set of the pulses. By alternating between the normal 
five-pulse primitives and the inverted hve pulses, cross¬ 
driving effects can be further reduced. (In fact, we note 


that this exactly eliminates all cross-driving that occurs 
via leakage in the VSM, because all pulses are always 
present at that distribution stage.) We refer to this as 
the symmetric 5-primitives technique and this is the tech¬ 
nique we implement in the main experiments described 
in Fig. 4 of the main text. Our simulations in Fig. |S10| 
show that the asymmetric 5-primitives technique already 
dramatically reduces the effect of cross-driving, and in 
the case of isolated single-qubit control, cross-driving is 
effectively eliminated completely using the symmetric 5- 
primitives scheme. This is also confirmed by measure¬ 
ments of cross-driving for the three selective-broadcasting 
schemes (Fig. S13), where we implement the symmet¬ 
ric 5-primitives technique. While we have only demon¬ 
strated this technique for the 5-primitives scheme, it 
could also be relatively straightforwardly applied to the 
sequential scheme, but the far better scaling of the 5- 
primitives scheme make it more interesting for scaling 
up to larger system sizes. 
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VII. CLIFFORD PULSE DECOMPOSITION 

The decompositions of the 24 single-qubit Clifford 
gates into a minimal set of 7r/2 and tt pulses and into 
the 5-primitives scheme are shown in Table |S3[ 


Clifford ID Minimal set decomposition 5-primitives decomposition 



First 

Second 

Third 

X ^/2 

to 

X ^/2 

x^ 

y-K 

1 

I 



0 

0 

0 

0 

0 

2 

U/2 

X ^/2 


0 

1 

1 

0 

0 

3 


Y .^/2 


1 

1 

0 

1 

0 

4 



0 

0 

0 

1 

0 

5 

Y-.n 

X -^,2 


0 

1 

1 

0 

1 

6 


Y .^,2 


1 

1 

0 

0 

1 

7 




0 

0 

0 

0 

1 

8 

Y .^/2 

X .^12 


0 

1 

1 

1 

1 

9 

X .„/2 

Y .12 


1 

1 

0 

0 

0 

10 




0 

0 

0 

1 

1 

11 

U/2 

X -^/2 


0 

1 

1 

1 

0 

12 

X -„,2 

U/2 


1 

1 

0 

1 

1 

13 

U/2 



0 

1 

0 

1 

0 

14 

X -„,2 



0 

0 

1 

1 

0 

15 

X ^ I 2 

to 

X -^,2 

1 

1 

1 

0 

1 

16 

Y-^/i 



0 

1 

0 

0 

1 

17 

X ^ I 2 



0 

0 

1 

0 

0 

18 

X ^/2 

U/2 

X ^/2 

1 

1 

1 

0 

0 

19 

Y .^,2 


0 

1 

0 

1 

1 

20 

X ^,2 

y^ 


1 

0 

0 

0 

1 

21 

X ^,2 

Y .^,2 

A 

to 

1 

1 

1 

1 

1 

22 

Ur/2 



0 

1 

0 

0 

0 

23 

X -^,2 

y^ 


1 

0 

0 

1 

1 

24 

X ^,2 

U/2 

1 

to 

1 

1 

1 

1 

0 


TABLE S3. Two decompositions of the 24 single-qubit Clif¬ 
ford gates. The first, taken from Ref. [S], minimizes the num¬ 
ber of 7r/2 and tt pulses around the ±x and ±t/ axes. The 
second is our decomposition into 5 primitives. Pulses are ap¬ 
plied from left to right. 


VIII. COMPILED SELECTIVE 
BROADCASTING ALGORITHM 


A. Finding the optimal pulse sequence 


When using a selective broadcasting architecture to 
send pulse sequences to multiple qubits, pulses can be 
directed to any subset of the qubits, but distinct pulses 
may not be applied simultaneously. In compiled selective 
broadcasting, the total number of pulses required to im¬ 
plement single-qubit gates on all qubits of a multi-qubit 
system is minimized by searching all possible combina¬ 
tions of single-qubit Clifford decompositions and group¬ 
ing together like pulses where possible. In this section, 
we introduce an algorithm for determining the short¬ 
est compiled pulse sequence implementing independent 
single-qubit Clifford gates on n qubits. 

On average, there are approximately 38 distinct de¬ 
compositions for each single-qubit Clifford gate, given the 
basis set of X and Y pulses: {J, Tjr, Aj-.,r/ 2 ) ^± 7 r/ 2 }j 
resulting in approximately 38" different decompositions 
for a given n-qubit combination of Cliffords. Here, we 
only consider sequences of up to four pulses, because the 
5-primitives decomposition identified in Table [S3| already 


provides a recipe for decomposing an arbitrary n-qubit 
Clifford combination into five pulses. We do not include 
trivial decompositions where sequential pulses cancel out. 

Given a particular choice of n Cliffords ,..., Ca ^), 
where Ui is the Clifford ID for qubit i, we write a spe¬ 
cific decomposition as ((Pi\ ..., J ,..., (Pf,..., P"^)), 
where Pj is the jth of rrii pulses which implement Cq,^. 
While this already fixes the order in which pulses must 
be applied to individual qubits, we still have the freedom 
to choose in which order the distinct pulses are applied to 
different qubits. For each possible decomposition, we use 
the following recursive algorithm to search and minimize 
over all possible pulse orderings. 

We first define an empty broadcasting sequence Pseq to 
store the compiled multi-qubit pulse sequence. In order 
to convert from parallel single-qubit pulse sequences to 
the single broadcasting sequence Pseq, we define a vector 
of indices (3 = (/3i,..., /3„) to store the current position 
in each single-qubit sequence. Initially, (3 = (1,...,1). 
At each instant, Pp = , ■ ■ ■, Pp ) contains the next 

pulses to be applied to each qubit. When Pi = rrii + 1, 
the pulse sequence for that qubit is completed, and so 
there is no Pp. to be added to Pp. The recursive part of 
the algorithm then proceeds as follows: 

1. Define Pp to be the set of distinct pulses in Pp. 

If:: Pp is empty, store the number of pulses in Pgeq 
and abort this recursion branch (Pgeq is a com¬ 
pleted pulse sequence such that all Cliffords 
are applied to the corresponding qubits). 

Else:: Continue. 

2. For each pulse P in P^, perform the following steps: 

(i) Append P to Pseq; 

(ii) Copy indices (3 to ...,/3))™); 

(iii) For all indices i for which Pp_ = P, increase 
the index Pf™ = /3j -|- 1; 

(iv) Recursively loop to step 1 using the new in¬ 
dices (3’^™ for (3. 

After considering all possible pulse sequences and loop¬ 
ing over all possible decompositions, choose the sequence 
with the minimum number of pulses Np in Pseq. 

This algorithm determines the minimum number of 
pulses Np required to implement a given n-qubit com¬ 
bination of single-qubit Clifford gates. However, this be¬ 
comes prohibitively resource intensive as the number of 
qubits increases. It is therefore important to assess how 
the performance of the optimal compiled decomposition 
compares with the 5-primitives decomposition, which can 
be applied to any number of qubits without any extra 
overhead in resources (in neither calculation time nor se¬ 
quence length). To do this, we use the average number of 
pulses (Np) required per n-qubit combination of Cliffords 
(per Clifford). 

In the case of compiled selective broadcasting, finding 
(Np) requires minimizing the sequence length for all 24" 
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Exact calculation 

Random sampling 

n 

(Np) 

n 

(Np) 

1 

1.875 

5 

4.139 (2) 

2 

2.925 

6 

4.380 (12) 

3 

3.521 

7 

4.570 (15) 

4 

3.874 

8 

4.721 (10) 

5 

4.137 

9 

4.808 (14) 



10 

4.857 (24) 


TABLE S4. The average number of pulses (A^p) required 
to perform one Clifford gate on each of n qubits in com¬ 
piled selective broadcasting. Exact values were obtained for 
1 < n < 5 in under two hours, using the outlined optimiza¬ 
tions. Approximate values were obtained for 5 < n < 10 
using random sampling. These results are plotted in Fig. 4 of 
main text. 

possible Clifford combinations. This problem again scales 
exponentially with n. For example, for n = 5 qubits, this 
requires 24® • 38® « 6.3 • 10^^ repetitions of the complete 
recursive search described above. Nevertheless, by em¬ 
ploying a number of optimizations desribed in the next 
section, we have exactly calculated (JVp) for 1 < n < 5 
qubits in under 2 hours. Using random sampling (find¬ 
ing the shortest pulse sequence for a random sample of 
Clifford combinations), we also approximated (JVp) for 
5 < n < 10. Exact and approximate results for n = 5 
agree. As shown in Table |S4[ the improvement offered 
by compiled selective broadcasting over the 5-primitives 
method is already less than one pulse per Clifford and 
continues to decrease rapidly. Considering how badly 
the resource overhead scales for increasing numbers of 
qubits for finding a compiled sequence, it is questionable 
whether compiling offers any significant benefit over us¬ 
ing the prescriptive 5-primitives approach when scaling 
up to larger system sizes. 


B. Optimizations for the Clifford compilation 
algorithm 

In the first optimization, we place an upper bound 
on the pulse sequence length. The upper bound is 
given by the minimum number of pulses found so far that 
can compile a given Clifford combination. At each stage, 
the algorithm checks if the sum of the pulses in Pgeq ^^nd 
all distinct pulses left is equal to or greater than 
If this is the case, a shorter combination of pulses using 
Pgeq is not possible, so we stop considering this sequence 
and proceed to the next one. Initially = 5, as the 
5-primitives method proves that there is always a decom¬ 
position of an arbitrary number of Cliffords into 5 pulses. 
Note that, as the limit moves down, the frequency 
at which the algorithm stops considering sequences in¬ 
creases. 

The second optimization relies on decompositions with 
fewer pulses being more likely to result in an optimal Clif¬ 
ford compilation. The decompositions of every Clifford 
are therefore arranged in ascending number of pulses. 
The first decompositions compared are then those with 
the minimum number of pulses; these have the highest 
probability of finding an optimal Clifford compilation. 
Even if an optimal Clifford compilation is not found, it 
is more likely that will be low. This optimization 
is especially effective in combination with the first opti¬ 
mization. 

The third optimization places a lower bound on 
the number of pulses. For a given Clifford combination 
(Caj,..., Ca„), is found by looking at the minimum 
number of pulses A^p previously found for all n—1 Clifford 
subsets. Since A^p for the n Cliffords can never be less 
than Np for any of the n — 1 Clifford subsets, the maxi¬ 
mum length of the n — 1 Clifford subsets therefore places 
a lower bound A^"^ on A^p for the n Cliffords. This means 
that if a pulse sequence is found whose length is equal to 
Af*^, it is an optimal Clifford compilation, and all further 
search is aborted. This is in contrast to the first opti¬ 
mization, where only the particular sequence of pulses is 
aborted upon reaching N'^^. Furthermore, as n increases, 
it becomes increasingly likely that the lower bound is 5. 
In this case, A^^'^ = and so the 5-primitives method 
is an optimal Clifford compilation. This optimization re¬ 
sults in the largest gain in computation time, by several 
orders of magnitude. 

In the fourth and most complicated optimization, all 
decompositions composed of three pulses or less are sepa¬ 
rated from those composed of four pulses. First, all com¬ 
binations of Clifford decompositions composed of three 
pulses or less are compared. This reduces the average 
number of decompositions per Clifford, from 38 to 7, re¬ 
sulting in an exponentially reduced number of total de¬ 
composition combinations. It is, however, not always 
the case that the optimal Clifford compilation is found 
using only up to three pulses per decomposition; some¬ 
times optimal Clifford compilations requires that one of 
the decompositions is composed of four pulses. However, 
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after comparing decompositions of three pulses or less, 
these four-pulse decompositions only need to be consid¬ 
ered when < 4 and = 5. If there is a sequence 
containing a four-pulse decomposition that outperforms 
any found using up to three-pulse decompositions and 
the 5-primitives method, the sequence must consist of 
four pulses. Only one Clifford then has a four-pulse de¬ 
composition, while all other Cliffords are subsets of these 
four pulses. We therefore loop, for every Clifford, over 
each of the four-pulse decompositions, and test whether 
every other Cliffords can be decomposed into a subset of 


these four pulses. This changes the comparison of four- 
pulse decompositions from scaling exponentially with n 
to scaling linearly. 

The fifth and final optimization is only of use when 
all different Clifford combinations need to be considered 
to determine (fVp). It stems from the observation that 
an optimal compilation for a certain Clifford combina¬ 
tion (Caj,..., is the same as for any permutation of 
those Cliffords. We therefore only determine an optimal 
Clifford compilation when /?i < • • • < /?«. This reduces 
the number of calculations exponentially (81 times fewer 
computations for n = 5). 
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